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Formulas for Inverse Osculatory Interpolation 

Herbert E. Salzer * 

% Formulas for inverse osculatory interpolation are obtained by inversion of Hermite's 
formula. They cover the cases for w = 2(l)7, where n is the number of points required in 
direct osculatory interpolation. The formulas provide an improved means for inverse inter- 
polation in the case where the first derivative is either tabulated alongside the function or 
is easily obtained. 



The author has previously given inverse interpo- 
lation formulas [1, 2] 2 for finding x=x -\-ph from 
f(x), in terms of f i =f(x i ), where the x t =Xo-\-ih are 
equally spaced at intervals of h. Those formulas 
were obtained by the inversion of Lagrange's inter- 
polation formula. As the Hermite osculatory inter- 
polation formula [3, 4, 5] in terms of the function 
and its first derivative at preassigned points x t is 
more accurate than the Lagrangian formula, inverse 
interpolation formulas obtained from the inversion 
of Hermite's formula would be expected to be more 
efficient than the corresponding formulas obtained 
from Lagrange's formula. 

There are many tables where the first derivative 
is either tabulated alongside the function (e. g., 
Bessel functions of the first or second kind, which 
give J ] (x)=-J^x) [6, 7] or Y 1 (x)= -Y^x) [8], 
respectively, or probability functions [9] ), or where 
the derivative is very easy to obtain (e. g., tables of 
the more elementary functions and their integrals, 
such as tables of sine, cosine, or exponential inte- 
grals [10, 11] ). In all such tables, and in many 



others, the user will find the inverse interpolation 
formulas given below to be particularly convenient, 
especially when the tabular interval h is too large for 
sufficiently accurate inverse interpolation, using the 
formulas in terms of the functional values alone. 

These present formulas give p in terms of f(x -\- 
ph)=f, f(xi)=ft, and f(x t )=fi for Xi=x +ih at 
equally spaced intervals h, where the i ranges from 
— [ (n-l)/2] to fn/2], for ?i=2(l)7, where n is the 
number of points required in direct osculatory 
interpolation. Although the direct interpolation 
formula for n=Q and n=7 is of the 11th- and 13th- 
degree accuracy, respectively, the inversion formula 
for p that is given below does not go beyond the 
10th- degree terms (and in most practical problems 
it is very rarely that one will go that far; in fact, 
the first few terms will usually suffice). Thus we 
do not use tin* coefficient of p n in the 6-point Hermite 
formula nor the coefficients of p l \ p l2 y and p ]3 in 
the 7-point Hermite formula. For every n, we 
define r=(f—f )/hf£, and corresponding_to n=2(l)7, 
quantities s, t, u, v, w, x, y, z, and z are defined 
as follows: 



n = 2 

^(- e SUo-fi}-h{2f +f 1 })/hf , 
t=(2Uo-fi}+h{fo+f 1 })/hf 0y 
u=v=w=x=y=z='z=0. 
n = 3 

«=(/.,-2/o+/i+| {/-i-/;))/V; 

<=(-f {/-!-/*}-£ {/-,+8/.+/iJ)/A/o, 

u={~\ {/- 1 -2/ +/ ] }-^ {/'_:-/;})/^, 
*=(! {/-i-/i}+i {/-.+4/;+/;})/¥;, 

w=x=y=z=z=0. 
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S= (l^8 {56/-i-297/ +216/ 1 +25/ 2 J+^ {2/!,- lS/J-lS/J-/;})/^, 

^=(~ {25C/_ 1 -/ 2 )+i35(y -/i)}+^ {/'- 1 +72/;+63/;+4/ 2 J) / /a/; ) 

»=(l58 {59/-,-54/ -27/ I + 22/ 2 }+A j 7/ '_ i+18/ ;_ 9/ ;_2^jy^ > 
V= (~lk {52/-i + 81/„-108y i -25/ 2 }-^ {5/'_, + 36/;+27/;+2/ 2 }) / //i/;, 



y=z=z=0. 

n=5 



s=(^[7(y_ 2 +/ 2 )^128(/_ 1 +y i )-270/ }+^{(/'_ 2 -/ 2 )+32(/'_ 1 -/0}) / /¥;, 

«=(-gJ i {3i(/- 2 -/ 2 )+i408(/_ I -/ I )j- I ^ i |(/'_ 2 +/ 2 )+64(/'_ 1 +/;)+36o/;}y/iy;, 

«=(-2^{41(/- 2 +/ 2 ) + 256(/_ I +/ 1 )-594/o}-^{3(/'_ 2 -/ 2 ) + 64C/'_ 1 -/;)})/^;, 
^=( TI ^{91(y- 2 -/ 2 ) + 1792(/_ 1 -/o} + r ^{3(/'_ 2 +^) + 128(/'_ 1 +/;)+396/;}yA/;, 

w= (ni { 13 (/^+^)+32(/- 1 +/ I )-9o/„J+A j(/'_ a -/ 2 )+i2(/'-i-y;)})/A/;, 

^=(-g^{29(/_ 2 -/ 2 ) + 272(/_ I -/ 1 )}-^{(/'- 2 +/ 2 ) + 24(/'_ 1 +/;) + 60/;}) / 'Vo, 

y=(-^{ii(/- s +/ 2 )+i6(/-i+/ 1 )-54/„J-^{(/L 2 -/D+8(/: 1 -/;)})^/; 

s=(^{25(/_ 2 -y 2 )+160(/_ 1 -yj}+Jg{(/'_ 2 +/ 2 )+16(/'-i+/0+36/i}yA/; 

2 = 0. 

n=6 
s=(^J76/_ 2 +2375/_ 1 -85OO/„+5000/ 1 +10O0/ 2 +49/ 3 } 

+ ^{3/'_ 2 + 150/'_ 1 -400^-600/;-75/ 2 -2/ 3 Jy^, 

t=(- ^^ f 3327/_ 2 -4- 1 69500/_ , - 8000/ - 1 68000/, + 2625/, + 548/i 

f 21/'_ 2 4-1500/'_ 1 +8600/^+1200/;-75/ 2 -4/ 3 Jy/t/;, 



3600 ^ 
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%=(-2j^oo{943/_ 3 -2050/_ 1 -62550/ u +45200/, + 17575f, + 882/ 3 } 

~~ 1416 { ] 3/l2 + 260/ - ' ~ 240 °/o- 296 °/i - 4 45/2 ~ 1 Vt \)/hf , 
»=(§6ioo } 5729 /-2+ 12140 °/-i- I6OOO/0- 1 1 8400/, + 6175/ 2 4- 1 096/,} 

+2^0 {37/'-,+1720/'-i+5140/5+ 1120/;-155/ 2 -8/ 8 }Va/ , 
w=( rj ^(l274/_ 2 -84875/_ 1 -189000/ +175000/ 1 +92750/ 2 +4851/ 3 } 

+14Im{ 36/ '- 2_2175/ '- 1-198 ^ 



r=(- Ii ^J5743/_ 2 +50125/_ 1 -22000/ -46000/ 1 4.10625/ 2 +1507/ 3 



1 39/'_,-l 925/:, + 1 900/;+ 100/;- 225/ 2 - 1 1/ 3 



4800 v 



y - (j^ J 1 72/_ 2 + 3275/_ , + 3900/o- 4600/, - 2600/, - 147, 

M80[' 
z=(gg^j351/_ 2 +750/_ 1 -4000/o+2625/ 2 +274/ a 



+ 4M I'"- 2 + 45/1 ' + 200/ ° + ] 80/; + 35/ 2 +/« 



+28^0 ] 3/1 ; + 30/ '" » ~~ 20/ "" T 2Q/ ' ; " W* ~ W* 



(- i 3^5J1066/_ 2 4-10625/_ 1 +9000/o-13000/ 1 -7250/ 2 -441/ a 



~ 14^00 1 8/ '- 2 + ] 75/1 ' + 600/ " + 500/; + ] 00/ 2 + 3 /a 
n = 7 
s= (3600o{ 157 C/-3+/ 3 )+4968(/_ 2 +/ 2 )4-43875(/_ 1 +/ 1 )-98000/ ( 

+ t4o{ c/ "- 3 ~"^ )+54 ^ 12 ^ )+675 ^ 1 '-/' ) 

*=(- T 0^oo{l670'-3-/3) + 8667(/_ 2 -/ 2 ) + 192375C/-i-/ 1 ) 

~«| c/ - 3+ - / ' )+81C/ '- 2+ - /2) + 2025 ^- i +/')+ 98() Q^ 



t= (~64W0"0 7339 C/"-3+/3) + 213786(/_ 2 +/ 2 ) + 631125Cf-,+/,)-1704500/ } 

"21600 { 47 (/'-3-/3) + 2403(/'_ 2 -/ 2 ) + 20925(/'_,-/;) 

^(i^owf 5206 ^- 8 ^ 3 ^ 249831 ^ 2 ^ 2 ^ 2517750 ^- 1 ^ 1 )} 

/) f 

+ J296Q0 94 (/'-3+/ 3 ) + 7209(/'_ 2 +/ 2 ) + 125550(/'_,+./;) + 340900/ 



w=(^ T g 1 t QQQJ52109(/_3+/3) + 1310016(/_ 2 +/ 2 )+1525875(/_,+/,)-5776000/ ( ,! 

A f J 

+ 172800 337 (/'- 3~/i) + 15648 (/'_ 2 -/ 2 ) 4- 88275 (/'_ , -/| ) 






)A/o, 



'A/S, 



)A/ , 



)/Vo, 

)/Vo, 
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= (~5l8iOOo{ 18493C/ - 3-/3) + 772368( - / - 2 ~- /2) + 4174125(/ - 1_ - /l) } 

^{337(/'_ 3 +/s)+23472(/'_ 2 +^)+264825(/: i +/;)+5776O0/ij)/A/; 



518400 



y=(- 172 g 000 J6217(/_3+y3) + 117408(/_ 2 +/ 2 ) + 72375(/_ 1 +/ 1 )-392000/ | 

-57feo{ 41 (/'-3-/3)+1584(/'- 2 -/2)+6675f/_ 1 -yD})/¥i, 

Vl728 



J 



728000 



2209(/_3-/ 3 ) + 70584(/_ 2 -/ 2 ) + 272625(/_ 1 -/ 1 ) 



+ 



172800 



^{4K/'_3+/;)+2376(/'_ 2 +/D + 20025(/'. 1 +/;)+39200/^jy/ l /; 



^=(^000 (2783 (/_ 3 +/ 3 ) + 38592 C/_ 2 +/ 2 ) + 14625 (/-i +/i) - 1 12000/oi 

+ i72^{ 19V '- 3 ~ /3)+576(/ '- 2 -- /2)+2025 ^'- 1 -- /:) })/^- 

For every n, the formula for p is the following: 

p=r-r 2 s+r*(2s 2 —t) + r*(-5s*+5st-u)+r 5 (Us*— 21s 2 t+3t 2 +6su-v)+r 6 (-42s 5 +84:sH 
-28st 2 -28s 2 u+7tu+7sv-w) + r 7 (l32s 6 -330sH+180s 2 t 2 +120s*u-12P^ 
+ 8^+86<w-z) + r 8 (-429s 7 +1287^-^ 

—4:5su 2 —90stv—4:5s 2 w+9uv+9tw+9sx—y)+r\US0s 8 ~5005sH+5005sH 2 +20Q2s 5 u 
-1430s¥-2860s 3 ^-715s 4 y+55£ 4 +660^+^ 

-110suv—110stw—55s 2 x+5v 2 +10uw+10tx+10sy—z)+r 1 \—4:8Q2s 9 +1944:8sH-240^ 
-8008s Q u+lOOlOsH d +15015sHu+3003s 5 v-l001s^-6006sH 2 u-2002s 3 u 2 -4:004sHv 
-l001s A w+286Pu+858stu 2 +858sHw+858st 2 v+858s 2 uv+ 
-66sv 2 —132suw—132stx-66s 2 y+nvw+llux+nty+llsz-z)+ • • .. 
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